Can the dynamics of an atomic glass-forming system be described as a continuous 

time random walk? 
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We show that the dynamics of supercooled liquids, analyzed from computer simulations of the 
binary mixture Lennard- Jones system, can be described in terms of a continuous time random walk 
(CTRW). The required discretization comes from mapping the dynamics on transitions between 
metabasins. This comparison involves verifying the conditions of the CTRW as well as a quantitative 
test of the predictions. In particular it is possible to express the wave vector-dependence of the 
relaxation time as well as the degree of non-exponentiality in terms of the first three moments of 
the waiting time distribution. 

PACS numbers: 



The dynamics of supercooled liquids is a very com- 
plex process with many non-trivial features such as non- 
exponential relaxation, decoupling of diffusion and relax- 
ation, significant correlated forward-backward processes 
(e.g. cage effect), and increasing length scales of re- 
laxation, just to mention some of the most prominent 
[H, 0, H[ ■ The complexity of the dynamics originates from 
the highly cooperative nature of the dynamical processes. 

Several phcnomcnological models have been proposed 
which attempt to describe the dynamics of supercooled 
liquids in relatively simple terms, thereby implying some 
kind of coarse-graining to get rid of the microscopic de- 
tails of the dynamics. In the free-energy energy 0, H, [E 0] 
and the RFOT models [8( the system relaxes, possibly 
in a multi-step process, between different states. One 
prominent example is the trap model [E @ , postulating 
a sequence of escape processes where the waiting time 
t in a configuration is fully governed by its energy and 
the new configuration is randomly chosen from the set of 
all possible configurations. Thus, the dynamics is fully 
described by the waiting time distribution <p(r). Ex- 
tending this model by the spatial aspects of the relax- 
ation processes one would, it its simplest version, end 
up with a continuous-time random walk (CTRW). Note 
that in general continuous-time random walks (CTRWs) 
[E [E EE, EH as well as the related Levy walks [H, EH 
are often used to describe anomalous dynamic proper- 
ties, characterized by non-trivial power-law behavior of 
quantities such as the mean-square displacement. 

In recent years the facilitated spin models have been 
revitalized to grasp the dynamics of supercooled liquids 
[II [IS EE 03 ■ They are thought to reflect the heteroge- 
neous mobility field of molecular glass- forming systems. 
One spin corresponds to a small volume which is either 
unjammed or jammed (spin up or down). The ability of 
a spin to flip is exclusively governed by the orientation 
of the adjacent spins. Self-diffusion has been introduced 
by postulating a random walk of the particle with the 
chance to move if the old as well as the new site is mo- 
bile [17J, UM> uM ■ This dynamics is also described in terms 
of a CTRW although for the model variant (East model) , 
supposed to describe fragile systems, a direct mapping is 



not possible 

Using the CTRW picture in the context of these phe- 
nomenological models does not necessarily imply that 
it is of relevance for microscopic glass-forming systems. 
Here we analyze a binary Lennard- Jones system (BMLJ), 
a standard model of supercooled liquids 20], via com- 
puter simulations. The goal of this work is threefold. 
First, we explicitly show that the required conditions for 
the applicability of the CTRW approach are fulfilled to 
a very good approximation. Second, we prove that the 
CTRW approach not only allows one to obtain the wave- 
vector dependent relaxation time (as already discussed, 
e.g., in [l8[) but also the non-exponentiality of relaxation. 
Third, the predictions are verified by explicit comparison 
with the numerical simulations. 

We analyze a BMLJ system with N=65 particles at 
T = 0.5 which is slightly above the mode-coupling tem- 
perature. It has been shown that this system is large 
enough to recover the diffusion constant without signif- 
icant finite size effect in the ran ge o f temperatures ac- 
cessible by computer simulatio ns I2lll22l| . Details of the 
model are described elsewhere 2oT 21 1. The discretiza- 
tion of the dynamics, required for the application of the 
CTRW approach, results from the use of inherent struc- 
tures, i.e. local minima of the potential energy land- 
scape [IE [24]], or the use of metabasins (MBs) [22T|. Of 
particular relevance in this work is the incoherent scat- 
tering function S(q,t) — (cos q[x(t + to) — x(to)]) where 
the brackets denote the average over all particles and 
all to- Furthermore, x(t) is the x-coordinate of a parti- 
cle. The first decay at short times to a value / < 1 is 
due to the fast /3-relaxation whereas the long-time relax- 
ation reflects the a-relaxation. It is often described by a 
KWW function / exp[-(t/T KWW )^ KWW ]. When analyz- 
ing S(q,t) for the sequence of inherent structures rather 
than actual configurations it turns out at temperatures 
close to the mode-coupling temperature that the short- 
time decay disappears and the decay is fully related to the 
a-relaxation with identical values tkww, Pkww [HI- 
Not surprisingly, the same holds for the sequence of MBs 
(data not shown). From now on, S(q,t) will represent 
the case of MBs, thereby describing the a-relaxation. 
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Figure 1: tti(x) for different values of the waiting times before 
the corresponding MB transition. Within statistical uncer- 
tainty tti(x) does not depend on the waiting time. 



Two important observables enter the CTRW approach: 
(i) the waiting time distribution <p{r), (ii) the probabil- 
ity 7Ti (x) that a particle during a transition between two 
MBs moves a specific distance along some fixed direction 
(here: x). More generally, ir n (x) expresses the corre- 
sponding probability after n MB transitions. The Fourier 
transform is denoted 7r„(g). Under conditions (C1)-(C3), 
which form the basis of the CTRW approach and are dis- 
cussed below, it is possible to express S(q, t) in terms of 
if (t) and TT]_(q). 

(CI) 7Ti (cc) does not depend on the waiting time since 
the previous transition. From the data in Fig. 1 the valid- 
ity of (CI) directly emerges. Only for the longest waiting 
times, which only have a very low probability (as reflected 
by the noise), minor deviations occur. As a consequence 
the spatial and temporal contributions separate to a very 
good approximation and one can write 

N 

S(q,t) = lim S N (q,t) = lim V] S n (t)ir n (q). (1) 

ri=0 

Here S n (t) denotes the probability to have exactly n tran- 
sitions during time t . This is the central equation of the 
CTRW because it expresses the total dynamics during 
time t in terms of discrete processes with well-defined 
probabilities. 

(C2) Successive waiting times are statistically uncor- 
rected so that the time evolution can be regarded as a 
sequence of randomly chosen waiting times. This has 
been already shown in Ref.[26|. Therefore S n (t) can be 
expressed in terms of the waiting time distribution tp(r) 
O d?! (see EqJS] below). Using the numerically deter- 
mined <p(t) and ir n {q), one can compare S(q max ,t), ob- 
tained from simulation, with the estimation Eq[T] where 
qmax is the maximum of the structure factor; see Fig. 2. 
The agreement is very good except for minor deviations 
for very long times. Of the order of 10 2 MB transition 
processes are required to have complete relaxation. 




Figure 2: Comparison of the actual incoherent scattering 
function S{q ma x,t) with the estimated function SN{q m ax,t) 
for different values of N. 

(C3) Subsequent transitions are spatially uncorrelated. 
The underlying Markov hypothesis can be formally writ- 
ten as 

7r n (x) = J dx'n n -i(x')ni(x - x'). (2) 

In Fourier-space this convolution reads 

n n (q)=ni(q) n . (3) 

From the analysis of the mean square displacement 
(MSD) in Ref. [2l| it became clear that there exist mi- 
nor forward-backward correlations for the MB transitions 
so that (C3) cannot hold in a strict sense. However, 
due to the expected locality of forward-backward tran- 
sitions one may expect that for longer length scales, i.e. 
smaller q, they become irrelevant. Indeed, one has a well- 
defined limit a 2 = lim„^ 00 (l/n)(a; 2 ) 7r == 0.005 which is 
slightly smaller than (x 2 ) 7ri — 0.009 [21| where (x 2 ) j 
j dxx 2 f{x). To check this in detail, we have analyzed 
the n-dependence of TT n (q), shown in Fig. 3 for different 
values of the wave-vector q. Interestingly, for q = q max 
the limiting behavior ir(q) = (T: n (q)) 1 ^ n — const is al- 
ready reached for n > 5, as reflected by the straight line. 
For smaller q- values Eq[3] holds even better. Since in the 
range of relevant q values one has a 2 q 2 <C 1 the term 
7r(g) can be approximated by 1 — q 2 a 2 /2. Using inher- 
ent structures rather than MBs the large-n regime would 
be only reached for n « 10 3 [2l| . This would strongly 
invalidate (C3). 

Using (C1)-(C3), and substituting all n n (q) by Tt(q), 
the temporal Laplace transform of the incoherent scatter- 
ing function, i.e. S(q,X), can be calculated analytically, 
yielding the Montroll- Weiss equation Q. Unfortunately, 
the inverse Laplace transform of S(q, A) cannot be an- 
alytically performed to calculate S(q,t). Therefore we 
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Figure 3: The n-dependence of n n (q) for different values of q. 



proceed in a somewhat different way. First, we define 



T (q)= / dtS(q,t) 



and 



J dttS(q,t)' 



(4) 



(5) 



To(q) denotes the relaxation time at wave vector q and 
(i m reflects the shape of S(q,t), based on the differ- 
ent moments. Whereas for exponential relaxation one 
has (3 m = 1 it decreases when S(q,t) decays in a non- 
exponential manner. In case of KWW relaxation one has 
(3 m = T 2 {l/p KWW )/{[3 KWW T{2/p KWW )) where T() de- 
notes the T-function (e.g. Pkww = 1/2 corresponds to 
An = 1/3). Pm depends monotonously on Pkww ■ Thus, 
[3 m is a measure of the degree of non-exponentiality. 

Our goal is to find simple expressions for To(q) and 
(3 m (q) . For this purpose one can introduce the persistence 
time distribution £(r). It denotes the probability that 
for a random starting point in time the next transition 
occurs a time r later [3, [H[ ■ It is related to the waiting 
time distribution via 



(0) 



£(r) = J dt' V {t')/{r) v . 
Furthermore, it is related to So(t) via 



So(t) 



dt'£(t'). 



(7) 



For n > the Laplace-transform of S n (t) is given by 

5 n (A) = C(A)V(A)"-V(T> ¥3 . (8) 

Straightforward calculation yields / dt S n (t) = S n (\ = 
0) = (t) v for n > 0. Note that for two functions, con- 
nected by f(t) — dt'g(t'), one obtains 



(t n )f 



(t n+1 ) s 
n + 1 



(9) 



This implies J dtSo(t) — (r)f, i.e. the average persis- 
tence time. Using again EqJS] it can be also expressed 
as (t 2 ) v / '{t) v . Note that (t)^/{t) v S> 1 for a broad 
waiting time distribution, reflecting large dynamic het- 
erogeneities. Equivalently, this means that the time un- 
til the first transition after a randomly chosen time takes 
much longer than the typical time {t) v between succes- 
sive jumps. 

Using EqQ] together with Eqs[3]and[9]one obtains [HI 



ro(q)/(r) v = jjf- + 



, 7Ti (g) 



{r)i 



(t) v l-7ri(g) (t) v 



q 2 a 2 



(10) 



We note in passing that (r)^ can be identified with the 
structural relaxation time r Q [18j . To determine the sim- 
ulated To(q) via integration over S(q,t) we have first fit- 
ted S(q, t) by a sum of two stretched exponentials and 
then performed the integration analytically. 

In Fig. 4a we show the comparison with the simulated 
data. We have used (t)^/(t) v = 27, as determined 
from the numerically determined waiting time distribu- 
tion. The g-dependence of Tp(q) is qualitatively similar 
to the data reported in [28| and [2i|. Note, however, that 
with the present definition of To(q) and the reference to 
the MB dynamics for the definition of (p(r) and a 2 a 
parameter-free prediction of the g-dependence becomes 
possible. At large q (q > q ma x) the system relaxes some- 
what faster because the effective value of a 2 increases due 
to the relevance of forward-backward correlations (see 
above). For smaller q, given by 1/q 2 ~ (r)f/ (t) v cl 2 /2, 
there is a crossover of To (q) from the q-independent large- 
q limit to the small-q limit r (q) — 2(r) v / '{q 2 : a 2 ) . Thus, 
for large dynamic heterogeneities, i.e. low temperatures, 
this crossover may happen at quite large distances [l8l |. 
Similarly, these non-trivial features are also reflected by 
a specific time evolution of the self-part of the van Hove 
function G s (x,t) [3(|. The deviations of G s (x, t) from 
simple diffusion have been analysed in detail in [311 ] . 

For the discussion of (3 m (q) we first rewrite Eqf5]as 



/J" 1 ( 9 ) - 1 = Tp(q)/T 2 (q) 



(11) 



with T/}(q) = J dttS(q,t) - r§(q). Following the 
standard derivation of the Montroll- Weiss equation one 
can show after a tedious but straightforward calcula- 
tion with the ingredients, presented in this work, that 
{d/dq)Tf}{q) = [27j], i.e. T p {q) = Tp. For the evaluation 
of Tp we choose the limit q — > 00 where S(q,t) = So(t). 
Following Eq[5]the first term equals (r 2 ){/2 whereas the 
second term is given by (t)?, i.e. 



T P (q) = ^ - <r)? 



(12) 



which directly reflects the width of the persistence time 
distribution. Note that via Eq[5] Tp involves the third 
moment of the waiting time distribution <p(t). Inter- 
estingly, the q-dependence of f3 m (q) is fully governed by 
r o (q) ■ Thus, the degree of non-exponentiality displays ex- 
actly the same crossover-behavior as the relaxation time. 
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Figure 4: The (/-dependence of (a) To(q) together with its 
estimation via Eq llUI and (b) of /3 m {q) together with its esti- 
mation. In the inset the validity of the theoretical expectation 
(d/dq)Tf3(q) = is tested. 



A comparison of EqsfTTI and [T2l with the numerical data 
is shown in Fig. 4b, showing again a good agreement. Ac- 
tually, due the extreme dependence of the third moment 
on the fine details of the long-time behavior of ip(r) a pre- 
cise estimation of Tp from tp(r) is not possible. Again, 
the deviations at large q reflect the more complicated dy- 
namics at short length-scales. The deviations for small q 
come from the trivial fact that Tp(q) results from a dif- 
ference of two very large numbers which, because of the 
nearly-exponential behavior, are very similar. 

In summary, the present work has shown that the 
CTRW approach, more or less explicitly used in different 
models of the glass transition, can indeed be numeri- 
cally derived for an atomic glass-forming system. This 
shows that after an appropriate coarse-graining proce- 
dure (here: the metabasins) the complex dynamics of 
supercooled liquids becomes relatively simple. Note that 
on a lower level of coarse-graining, namely the inherent 
structures, (C3) and thus the CTRW approach is strongly 
violated. In contrast, more coarse-graining, e.g. by join- 
ing some successive MBs, would start to change S(q, t) by 
rendering it more exponential. In analogy to the previous 
model considerations the CTRW approach is formulated 
for a subsystem of a large macroscopic system (one coop- 
eratively rearranging region, one probe molecule, here: a 
small system with periodic boundary conditions). Gen- 
eralization to large systems, thereby keeping information 
about possible correlations and predicting multi-point 
correlation functions, is a challenge for the future. 
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